Wave Scattering through Classically Chaotic Cavities in the Presence of Absorption: 

An Information-Theoretic Model 
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We propose an information-theoretic model for the transport of waves through a chaotic cavity 
in the presence of absorption. The entropy of the S-matrix statistical distribution is maximized, 
with the constraint ^TrS'S'^) — an: n is the dimensionality of S, and < a < 1, a = 0(1) meaning 
complete (no) absorption. For strong absorption our result agrees with a number of analytical 
calculations already given in the literature. In that limit, the distribution of the individual (angular) 
transmission and reflection coefficients becomes exponential -Rayleigh statistics- even for n = 1. 
For n 3> 1 Rayleigh statistics is attained even with no absorption; here we extend the study to 
a < 1. The model is compared with random-matrix-theory numerical simulations: it describes the 
problem very well for strong absorption, but fails for moderate and weak absorptions. Thus, in the 
latter regime, some important physical constraint is missing in the construction of the model. 



PACS numbers: Q5.45.-Hb,42.25.Bs,41.20.Jb 



Systems involving multiple elastic scattering of any 
kind of waves (sound, microwaves or light) -whose in- 
terference leads to strong fluctuations in the transmitted 
intensity- can be described in very much the same way 
as electron systems. However, there exists an important 
difference: the interference pattern for classical waves 
can be affected as a result of loss or absorption, which 
is absent in electron systems. The issue is extremely im- 
portant from an experimental point of view, because ab- 
sorption is always present and is often very strong. For 
diffusive transport the problem was intensively studied 
both experimentally and theoretically (|-^. The is- 
sue has also attracted attention in connection with the 
phase-coherent reflection of light by a disordered medium 
which amplifies radiation and the study of the rela- 
tion between absorption and dephasing ^,0. 

The analytical evaluation of the reflection-matrix sta- 
tistical distribution for a semi-inflnite disordered waveg- 
uide was performed, for arbitrary absorption, in Ref. M, 
and that of the S'-matrix distribution for a chaotic cav- 
ity with absorption and one propagating mode in each of 
two waveguides, in Ref. 0|. 

In the present paper we take up again the problem of 
the propagation of scalar waves traveling through a cav- 
ity -whose classical dynamics would be chaotic- in the 
presence of absorption, connected to the outside through 
a number of waveguides supporting an arbitrary number 
of propagating modes. Motivated by the success of an 
information-theoretic approach ||] to the study of chaotic 
scattering through cavities, we propose below an exten- 
sion of such models to study the effect of absorption. 
Since these models are based on the idea of doing statis- 
tics directly on the S matrix of the system -on the basis of 
the information which is physically relevant for the prob- 



lem in question- we believe that the present approach 
complements the analytical derivations mentioned above. 
We show that the two approaches agree in the limit of 
strong absorption, while for moderate and weak absorp- 
tion some relevant information is missed in our model. 

The scattering of waves through a cavity can be de- 
scribed by an S matrix that relates incoming and out- 
going amplitudes. The dimensionality n of the matrix is 
the total number of channels in all the waveguides. For 
two A'^-channel waveguides, n = 2N and the S matrix 
has the structure 



r t' 
t r' 



(1) 



where r, r' and t, t' are the A^-dimensional matrices of re- 
flection and transmission amplitudes with incidence from 
either waveguide. 

In quantum-mechanics, the universality classes for S 
matrices were introduced by Dyson . In the absence 
of any symmetries, the only restriction on S is unitarity, 
SS^ — I, due to flux conservation (the unitary or /3 = 2 
case). In the orthogonal case (/? = 1), S* is symmet- 
ric because of either time-reversal invariance (TRI) and 
integral spin, or TRI, half-integral spin and rotational 
symmetry. In the symplectic case (/? = 4), S" is self-dual 
because of TRI with half-integral spin and no rotational 
symmetry. The intuitive idea of equal- a-priori probabil- 
ities is expressed mathematically by the invariant mea- 
sure on the matrix space under the symmetry operation 
for the class in question, giving the circular orthogonal, 
unitary and symplectic ensembles (COE, CUE, CSE). 

Consider the orthogonal case. The potential appear- 
ing in the Schrodinger equation is real, with a strength 
uq, say. The resulting S matrix is unitary, SS^ — I, 



1 



and symmetric, S = S"^ . Suppose we analytically con- 
tinue Mo to complex values: uq = u'q — iu'^, the sign of 
the imaginary part ensuring absorption. Because of loss 
of flux, the resulting S is now subunitary, meaning that 
the eigenvalues of the Hermitian matrix h — SS^ lie in 
the interval between zero and one. However, the sym- 
metry property, S = S'^ , is not altered. We still speak 
of the orthogonal {(3=1) case. For the unitary one, a 
similar analytic continuation gives a subunitary and in 
general non-symmetric S matrix. In the scattering prob- 
lem of scalar classical waves, the orthogonal case is the 
physically relevant one. However, we deal below with 
both (3=1 and 2, the unitary case being presented as a 
reference problem, as it is often simpler to treat mathe- 
matically than the orthogonal one. 

Following Rcf. (pp. 63, 64), we introduce a uni- 
form weight in the space of sub-unitary matrices as 

= C0{I - SS^) H dxatdyab, (2) 

where Sab — Xab + iUab] Ha b '^^^'^ elements in the 
unitary case, but, in the orthogonal one, only over a < b. 
In this equation, the step function 9{H) (for a Hermitian 
H) is nonzero for iJ > (i.e. for H positive definite, 
so that all its eigenvalues are positive definite) and thus 
selects sub-unitary matrices. 

A complex n x n sub-unitary matrix can be written in 
the polar representation (Refs. pp. 63,64 and 0) as 

S = UDV. (3) 

The unitary matrices U and V are arbitrary in the uni- 
tary case, while V = U'^ in the orthogonal one. The 
matrix D is diagonal, with the structure Dab — Pa^ab 
(a = 1 , • • • , n) , with < pa < 1 • 

The explicit expression of the above measure (^) in 
terms of the independent parameters of the polar repre- 
sentation (|) is (Rcf. 0, pp. 63, 64) 

n 

a<b c 

where dfi{U), dfj,{V) are the invariant measure for the 
unitary group in n dimensions; dfi{V) is absent in the 
orthogonal case (3=1. 

More general statistical distributions of sub-unitary 
matrices carrying more information than the equal-a- 
priori probability (|^) can now be constructed using 
dfJ.subiS) as a starting point and writing 

dP{S)=p{S)dnsubiS). (5) 

In what follows we propose an information-theoretic cri- 
terion to choose p{S). The information-theoretic en- 
tropy S of the 5'-matrix distribution jl^ , S[p{S)] = 
— J p{S) \Tip{S)d^subiS), is maximized subject to the con- 
straint of a given average strength of the absorption. 



Mathematically, this is expressed by the average depar- 
ture from unitarity of our S matrices; we thus write the 
constraint as 

(TrSS^) =an, < a < 1. (6) 

Thus a = corresponds to complete absorption and 
a = 1 to lack of absorption. We find 

dpi^\S)=Ce~^^^-'^'df,Us), (7) 

where the constant C and the Lagrange multiplier v en- 
sure normalization and the fulfillment of the constraint 
. The limit of no absorption, a = 1, is attained when 
the Lagrange multiplier v ^ —oo and the distribution 
concentrates on the unitarity sphere. The limit of com- 
plete absorption, a — > 0, is attained when v +oo: the 
distribution then becomes a (5-function at the origin and 
there is no exit signal. The result of Eq. (|^) -Laguerre 
ensemble for the variables p^- coincides, for strong ab- 
sorption, with that obtained in Ref. |^ for the diffusive 
waveguide and in Ref. ||^ for cavities with N = 1. 

The "ansatz" (|^) entails a number of properties and 
restrictions. ^From (Q), (||) and the properties of d^{U) 
|l3| , we see that, under the distribution (^, the average 
{S) = 0. Therefore, applications of the model (^ should 
be restricted to cases where prompt processes are absent 
and so (S) = |ll|,|l|,§. We shaU see at the end of the 
paper that for uniform volume absorption (in quantum 
mechanics, a constant potential —iW throughout the 
cavity) the S matrix can be obtained from that without 
absorption, evaluating it at the complex energy E -\- iW, 
i.e. S{E + iW). The so-called analyticity-ergodicity re- 
quirements (^S^li,^ ■ ■ ■ S'";^^) = {Sa^b^^' ■ ■ ■ {Sa^bS" 

are discussed in Ref. for unitary matrices: the same 
argument applies here as well, and the necessity to fulfill 
them follows. In the present case, their fulfillment follows 
from (|), (|), (|) and the properties of d^i{U) |l3|. 

As a brief excursion into the situation when {S) ^ 0, 
we note that, as W increases, S{E -t- iW) (S), which, 
in turn, is the so called "optical S matrix", that signals 
the presence of prompt or direct processes |l^Jl^,^. As 
a consequence, S, for large absorption, does not tend 
to zero, but to (S). With no absorption, the spHtting 
S = (S) -\- S^^ describes the problem in terms of two 
responses, associated with two distinct time scales: the 
prompt and the equilibrated one. It is then expected 
that strong absorption will affect the former response, 
that corresponds to short trajectories, much less then 
the latter one, that arises from long chaotic trajectories. 

We analyze some of the consequences of the ansatz 
(0) and, at the end, compare them with the results of 
random-matrix-theory (RMT) numerical simulations. 

a. The n = 1 case. This case, which describes a cav- 
ity with one waveguide supporting only one open chan- 
nel {S is thus the reflection amplitude back to the only 
channel we have), is, within our model (^, independent 
of the universality class /3. Eq. (||) for S in the polar 
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representation reduces to S = pexpit/; p 
reflection coefficient R. The uniform weight ^ 
distribution (Q) reduce to 

dpsubiS) = pdpde, dP{S) = Ce-'p'pdpde. 
The i?-probability density is 



^ represents the 



and the 



(8) 



D and v being given by 



< i? < 1, 



1 



(9) 



(10) 



For weak absorption, a sa 1, — > — oo and the distri- 
bution (^) becomes strongly peaked around i? = 1, i.e. 
the unitarity circle, reducing to the one-sided delta func- 
tion i5(l — i?) as a ^ 1. In the other extreme of strong 
absorption, v — > -|-oo, a k, D k, \/a and 



w{R) « {R) 



-1 ^-R/W 



(11) 



Rayleigh's distribution, with the average (i?) = a. 

b. The orthogonal case for arbitrary n. Using the re- 
sults of Ref . we find, for the average of an individual 
(angular) transmission or reflection coefficient 

{Tabf^ - {Ra^bf^ - al{n + 1) = (1/2) {R^af^ (12) 

We see the occurrence of the familiar backward enhance- 
ment factor 2. For the second moments we find 



(1) 



E 

a=l 



{\u,^\'\u,^\')M) 



(1) 



+2 



(13) 



The indices 1 and 2 indicate any pair of different chan- 
nels; (• • ■)q stands for an average with respect to the in- 
variant measure of the unitary group ||l^. For (-Raa)'^' 
one sets the two indices 1 and 2 equal. 

For a two-waveguide problem with one channel in each 
waveguide the S matrix is two dimensional [n = 2N = 

2). Eq. (H) gives (T)'^' = a/3, (R)^-^^ = 2a/3. Restrict- 
ing ourselves to the limit of strong absorption, a <|; 1, the 
Lagrange multiplier i/ = 3/2q; and we obtain 



(1) 



(T) 



(1) 



(14) 



Although we have not calculated the statistical distri- 
bution of T, result (^4|) is consistent with the Rayleigh 
distribution for the transmission coefficient T. 

For a large number of channels, n = 2N ^ 1, Eqs. 
([|) and (0) give Q 



{T!bf' = {RUbr'-^2^^2 



(Tab)^'^ 



(15) 



and similarly, (-R^^) 



(1) 



{Raa) 



(1) 



the relation be- 



tween first and second moments for an exponential dis- 
tribution with the average value ([l^), which becomes 
smaller as the absorption increases. 

For no absorption one reaches, in the limit n ^ 1, a 
Rayleigh distribution for Raa- Rcf. jl^ shows that, for 
the COE (« being an approximation for large n) 



w 



(1) 



{Raa) = C {1 — Ra 



(Raa) 



c. The unitary case for arbitrary n. In the unitary 
case, the statistical properties of a transmission coeffi- 
cient is identical to that of a diagonal or off-diagonal 
reflection coefficient. The following equations are thus 
written for Tab- We find, for its average, {Tab) = a/n; 



we have used the result |13| ( \Uia 



^ 1/n. The dif- 
ference in expectation values between the two symmetry 
classes is thus 



{Tabf^ - {Tab) 



(2) 



-a/[n{n+l)]. 



(16) 



For the second moment of Tab we have 

n2 \(2) 



a-y 

= [2/{n + 1)2] {[{n l)/n] {pjplf^ + {2/n) {p\f^] . (17) 



We have used the resufi ^ {\Uic\'^\Ui 
{1 + Saj) / [n{n + 1)]. 

For a two-waveguide problem with one channel in each 
waveguide {n — 2N = 2) we have (T) = a/2. In the 
limit of strong absorption, a <^ 1, the Lagrange multi- 
plier = 2/a and we obtain a relation like ( |l4|). For a 
large number of channels n — 2N ^ 1, Eq. (|17|) gives 
a similar relation, now for Tab, which is again consistent 
with Rayleigh's distribution. 

For CUE (i.e. no absorption) one reaches a Rayleigh 
distribution for n ^ 1. Ref. [|5| finds for the the distri- 
bution of a single transmission coefficient as 

w{Tab) ^C{l-Tabj"-' « (r„fc)-ie-^-/<^->. (18) 

d. Comparison with numerical simulations Some 
of our predictions are compared below with RMT nu- 
merical simulations for (3 = 1. The S matrices are 
constructed as S{E) = - [I,, - iK{E)]~^ [1^ + iK{E)], 

with Kab{E) = Y.xl>-anXb/ [E\ ~ E] and {In)ab = 5ab 

{a,b = 1,- ■ ■,n). The E\s are generated from an 
"unfolded" zero-centered GOE with average spac- 
ing A. The 7Aa's are statistically independent, real, 
zero-centered Gaussian random variables. At i? = 0, 

{Sab) = -[l + 7r(7L)/A]"'[l-7r(7L)/A]'^afc, and 
we require {S) = 0. In the quantum case, addition of 
a constant imaginary potential —iW inside the cavity 
makes the E\s complex and equal to Ex — iW (see also 
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Ref. 0). This is equivalent to evaluating the above ex- 
pressions at the complex energy E + iW, which makes 
S{E + iW) subunitary. Although Eq. (|^) gives the prob- 
ability distribution for the full S matrix and arbitrary n, 
we only analyze below individual (angular) reflection and 
transmission coefficients, for n — 1,2. 




FIG. 



R R 
1. The distribution of the reflection coefficient R for a 



cavity with one waveguide, supporting one channel: (left) in 
the presence of strong absorption {{R) = a = 0.034); (right) 
for moderate absorption ((7?) = a — 0.410). 
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FIG. 2. The distribution of the transmission coefficient T 
for a cavity with two waveguides, each supporting one open 
channel, in the presence of strong absorption (a — .049). 



Fig. 1 shows for n = I the results of the RMT nu- 
merical simulations (as histograms), compared with the 
present model for the corresponding value of a (continu- 
ous curves) . For strong absorption the model works very 
well, the agreement with Rayleigh's law being excellent, 
while for moderate and weak absorptions the model fails. 

Fig. 2 shows the distribution of the transmission co- 
efficient T obtained from a RMT simulation for n = 2. 
The agreement with the Rayleigh distribution with the 
centroid (T) = a/3 is excellent; (R) was also checked and 
found to agree with 2a/3. 

That individual transmission and reflection coefficients 
attain a Rayleigh distribution for strong absorption can 
be understood as follows. Sab{E + iW) coincides [pjT2| 
with the energy average of SabiE) evaluated with a 
Lorentzian weighting function of half-width W. If F'^"'"'" 
is the correlation energy, W can be thought of as contain- 
ing ^ m — W/V^"^^ independent intervals. If m 1, by 
the central-limit theorem the real and imaginary parts 
of Sab attain a Gaussian distribution, and |5afc| an ex- 
ponential distribution. This seems to be the situation 
captured by the maximum-entropy approach. 

Summarizing, the results presented in this paper in- 
dicate that wave scattering through classically chaotic 
cavities in the presence of strong absorption can be de- 
scribed in terms of an information-theoretic model. 

We benefited from the constructive criticism of the first 
version of the paper by C. W. J. Beenakker and P. W. 
Brouwer. One of the authors (PAM) acknowledges par- 
tial financial support from CONACyT, Mexico, through 
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